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d • Fermi coordinates are directly constructed in de Sitter and Godel spacetimes and 



the corresponding exact coordinate transformations are given explicitly. The quasi- 
inertial Fermi coordinates are then employed to discuss the dynamics of a free test 



' particle in these spacetimes and the results are compared to the corresponding gen- 



eralized Jacobi equations that contain only the lowest-order tidal terms. The do- 
main of validity of the generalized Jacobi equation is thus examined in these cases. 



1/-^ ' Furthermore, the difficulty of constructing explicit Fermi coordinates in black-hole 

o ■ 

■ spacetimes is demonstrated. 

7'. 

tl' PACS numbers: 04.20. Cv 



I. INTRODUCTION 

To interpret measurements in a gravitational field, access to locally inertial coordinates 
is indispensable. At each event in spacetime, the spacetime is locally flat; therefore, it is 
possible to introduce Riemann normal coordinates that constitute a geodesic system of co- 
ordinates that is inertial at the event under consideration. The Riemann normal coordinates 
are in general admissible only in a certain region around the event such that every point in 
this region can be connected to the event by a unique geodesic. 

In general, the physical interpretation of measurements by a free observer necessitates a 
continuous locally inertial system along the worldline of the observer. This can be achieved 
by Fermi coordinates, namely, a normal geodesic coordinate system in a cylindrical region 
about the worldline of the observer. Fermi normal coordinates 111 are the natural extension 




FIG. 1: Schematic construction of Fermi coordinates at P. 

of Riemann normal coordinates and play a basic role in general relativity theory P]. The 
notion of Fermi coordinates was first implicitly introduced in Ref . jl| ; a recent commentary 
on Fermi's pioneering work is contained in Ref. P|. 

Imagine a spacetime region with coordinates = [t, x*) and a reference observer O 
following a worldline (t, x*). We use gravitational units such that c = G = 1 throughout; 
moreover, we assume that the spacetime metric has signature +2. The observer carries an 
orthonormal parallel-propagated tetrad frame A'^^.^-j along its path such that A'^^^q-j = dx^/dr, 
where r is the proper time along the worldline of O. Thus A'^^q-j is the timelike unit vector that 
is tangent to the worldline of O and acts as its local temporal axis. Moreover, A^^^.^, i = 1, 2, 3, 
are orthogonal spacelike unit gyro axes that form the local spatial frame of the observer. 
At each event Q{t) along the worldline, consider the class of spacelike geodesies orthogonal 
to the worldline; these form a local hypersurface. Let P be an event with coordinates x'^ 
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on this hypersurface and consider the unique spacehke geodesic segment from Q to P. We 
define the Fermi coordinates of P to be = (T, X*), where 

T = T, X' = aeX^^'\ (1) 

Here is the unit vector tangent to the spacehke geodesic segment at Q and a is the proper 
length of this segment from Q to P as in Figure 1. Thus the reference observer O is always 
at the spatial origin of Fermi coordinates. 

The Fermi coordinates are in general admissible in the sense of Lichnerowicz jj] in a finite 
cylindrical region about the worldline of O. The spacetime metric in Fermi coordinates is 
given by 

goo = -l-^Ro^OJ{T)X'X^ + ■■■ , (2) 
go^ = -l''RoJ^k{T)X^X'' + ■■■ , (3) 
g^J = % - ^^R^kAT)X'X' + ■■■ , (4) 

where 

^RaPysiT) = _R^i,p^A'^(^)A''(^)A^(^jA'^(5) (5) 

is the projection of the Riemann tensor on the observer's tetrad along the reference tra- 
jectory. The Taylor series in Eqs. (0)^(111) can be expressed as g^^, = rj^y + h^y{X), where 
h^j,u{X) is a perturbation that can be expressed as a series expansion in powers of the spatial 
distance away from the reference trajectory. The nature of the infinite series in Eqs. ((2}- 
(j3)) has been discussed by a number of authors. The second-order terms given explicitly 



in Eqs. (j2I)^0 were worked out in Refs. ^ , la] . The third-order terms were first given in 

Ref . ItII and the fourth-order terms in Ref . S I ; moreover, these higher-order terms have been 

fl 

recently discussed in Ref. j9|. In the weak-field limit, the infinite series in Eqs. have 
been studied in Ref. 1^ . 

Consider the motion of a free particle in the Fermi coordinate system. The geodesic 
equation of motion 

__ + Y^^^uW = 0, (6) 

where If^ = dX'^/ds, can be written in terms of the modified Lorentz factor F = dT /ds and 
V = dX/dT based on the decomposition = r(l, V). On the other hand, the timelike 
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condition U^Uu = —1 implies that 



-^ = gm + '^g^^V' + 9^JV'VK (7) 
The equation for V(T) is the reduced geodesic equation given by 

^ + (r„,-r,,^)^^-o. (8) 

It is useful to consider the class of static observers in the Fermi coordinate system. These 
observers are generally accelerated with four- velocity Ug = {—goo)^^^^o in Fermi coordi- 
nates. The Lorentz factor of the test particle with four-velocity = r(l, V) with respect 
to the static observers is 

V^o-^)r. (9) 



In this paper, we are interested in the comparison between the consequences of these 
equations in the case of exact Fermi coordinates — given explicitly in de Sitter and Godel 
spacetimes — with equations based on the lowest-order tidal terms given explicitly in 
Eqs. ©-(HI)- In the latter case, Eq. (0) implies that 



(10) 



and the reduced equation of motion is the generalized Jacobi equation 
(PX' 



+ ^{'i'Rok.oV'V'' + ^R^kJlV^V' + ^i?ofc,7^W)X^' = 0. 
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(11) 



This equation is expected to be valid for |X| sufficiently small compared to a certain radius 
of curvature of spacetime. We assume the following initial conditions for Eqs. (jH)) and 
throughout this paper: At T = 0, X = and V = Vq such that |Vo| < 1. Let us note that 
for T > and away from the origin of Fermi coordinates, |V| could then in principle exceed 
unity along a timelike geodesic worldline. 

In some situations — for instance, along the axis of rotational symmetry of a Kerr black 
hole — one- dimensional motion is allowed. In general, the symmetries of the Riemann tensor 



5 



imply that for motion along the Z direction, say, Eq. reduces to 

^ + «:(1-2Z^)Z = 0, (12) 

where k(T) = ^Rtztz(T) and Z = dZ/dT. This equation approximates the reduced geodesic 
equation for \Z\ sufficiently small compared to |k|^^/^. In Eq. (fT^ . the critical speed of 
l/-\/2 ^ 0.7 should be noted: Motion at the critical speed is uniform, i.e. Z = constant ± 
T/v^, and the character of the motion changes depending on whether the initial relative 
speed is above or below the critical speed. The critical speed in gravitational motion is 
briefly reviewed in Appendix A. 

In a recent series of papers, we have investigated the general equation of motion of a 



particle in Fermi coordinates using the generalized Jacobi equation The results are 

particularly interesting if the speed of the test particle relative to the reference particle 
exceeds the critical speed given by \4 := 1/a/2. The astrophysical implications of these 
results have been worked out in Refs. jl^ . [ This general approach has been extended 
to the motion of charged particles in Ref. 1^. A major shortcoming of these studies is 
that only the first few terms of the series in Eqs. have been taken into account. It 

is therefore important to know the domain of validity of these astrophysically significant 
results. In fact, as explicitly demonstrated in Appendix B, the construction of the exact 
Fermi coordinates in black-hole spacetimes is a daunting task. Moreover, the transformation 



13|. For recent 



H-> remains implicit as well as approximate in these studies [11|, ll^ 
efforts to construct such approximate transformations explicitly, see Ref. 

To avoid the difficulties encountered in black- hole spacetimes ( cf. Appendix B), we turn 
to spacetimes with more symmetries. It is clear that explicit Fermi coordinates can only be 
constructed in rather special circumstances; however, the generalized Jacobi equation can be 
employed in general. What is the extent of agreement between these approaches? To answer 
this question, we choose two spatially homogeneous spacetimes: de Sitter's spacetime with 
ten Killing vector fields and Godel's spacetime with five Killing vector fields. 

The main purpose of the present paper is to study the infinite series in Eqs. ©^Q and 
the exact transformation ^ to Fermi coordinates explicitly in de Sitter and Godel 
spacetimes. The motion of a test particle is then studied using the exact Fermi metric and 
the results are compared to the implications of the generalized Jacobi equation (fT^ . where 
in de Sitter (Godel) spacetime k is a negative (positive) constant. Eq. can be integrated 
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in these cases; the behavior of the solutions have been studied in Ref. [hI. The comparison 
of the exact geodesic equation of motion in Fermi coordinates with the generahzed Jacobi 
equation makes it in principle possible to determine the extent of validity of the latter 
equation in these cases; moreover, a general mathematical treatment of such a comparison 
is presented in Appendix C. 

II. FERMI COORDINATES IN DE SITTER SPACETIME 

Consider de Sitter's metric in the (inflationary-model) form 

ds'^ = -df + A\t)6,jdx'dx^, (13) 



where A{t) is given by 



A{t) = e^* (14) 



and if > is a ("Hubble") constant. This metric satisfies Einstein's matter-free equations 
such that = A/3, where A > is the cosmological constant. The spacetime region 
x'^ = (t, x) described by equation ((T)) is only part of de Sitter's spacetime of constant 
positive curvature: a detailed historical description of spacetimes of constant curvature is 
contained in Ref. 

The geodesies of the metric (fTSj) are obtained from the extrema of J ds. It is a simple 
consequence of this fact that A^dx^ds, i = 1,2,3, are constants along the geodesies. Let 
us imagine a class of reference observers that are fixed in space, i.e. they have constant 
spatial coordinates x*, i = 1,2, 3. It turns out that these are free observers each with an 
orthonormal tetrad frame A'', ^ such that 

(a) 

where t = r is the proper time along the worldline of the observer. Moreover, it is simple 
to verify that this tetrad frame is parallel transported along the geodesic worldline of the 
observer. 

The projection of the Riemann tensor on the tetrad of a reference observer can be 
expressed as a 6 x 6 matrix 71 = {71 ab), where A and B are elements of the set 
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{01, 02, 03, 23, 31, 12}. We find that 



7^ 





-/ 



(16) 



where / is the 3x3 unit matrix. 

We need to find the general solution for spacelike geodesies. To this end, the spacelike 
geodesic equation in this case reduces to 

dx^ a ^^^^ 

(18) 



da A^{t)' 
dtV 



- 



da J A^{t) 

where a is the proper length along a spacehke geodesic, Q, i = 1,2, 3, are constants and 



(19) 



The general solution of Eq. (|T8|) is 



e^* = Ccos{Ha + ri), 



(20) 



where rj is a constant. It follows from Eq. ()20|1 that Eq. ()17|1 can be solved exactly and the 
result is 



X 



tan{Ha + 7]) + Di 



(21) 



where Di, i = 1,2, 3, are constants of integration. 

To establish a Fermi coordinate system, we need to choose a specific reference trajectory. 
For the sake of simplicity and with no loss in generality, we choose the reference observer 
O : (t, x) = (t, 0). In the construction of Fermi coordinates (see Figure 1), the spacelike 
geodesic segment from Q to P is given by equations and (PT|) such that cr = at 
Therefore, 



e^'^ = C cos T], 



tanr/ + Di = 0. 



(22) 



Moreover, C,^ is the vector tangent to the spacelike geodesic segment at Q and is given in 
(if:, x) coordinates by 



e 



tani] 



' cos^ rj 



(23) 
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It follows from ^/^A'^j^q-j = that tanrj = and hence Di = 0, i = 1, 2, 3, from Eq. (j22I)- We 
choose T) = for simplicity; hence, C = exp(ifr). Finally, we compute the spatial Fermi 
coordinates of P using Eq. (^. It follows that event P has Fermi coordinates 

T = r, X = ae-^^C. (24) 
It is useful to recognize that a is the radial Fermi coordinate R, 

R = y^6~J0Xl, (25) 
so that we find from Eqs. (j^UI) and (PT|) that the transformation {t, x) i-^ (T, X) is given by 

e^' = e^^ cos{HR), (26) 

x = e-^^^^^^X. (27) 
HR 

De Sitter's metric in Fermi coordinates is ds"^ = g^^dX^dX^ , where 

goo = -cos\HR), (?o* = 0, (28) 

The Fermi coordinates cover a static spacetime region and are admissible for < i? < 
Tx/{2H)- in fact, ^oo = for HR = 7r/2. 

The form of the metric ()28|l - p9|l can be simplified if we introduce Fermi polar coordinates 

via 

X^ = R sin 9 cos 0, X^ = R sin 9 sin 0, X^ = R cos 9. (30) 
The Fermi metric in these spherical polar coordinates is given by 

ds^ = -cos\HR)dT^ + dR^ + ^sm\HR){d9^ + sin^ 9 d(j)^). (31) 

H 

It is interesting to note that this form of de Sitter's metric already appeared in de Sitter's 
original investigations, namely, in his 1917 paper on the curvature of space (cf. Eq. (2.5) of 
Ref. j3| and the discussion therein). 

Writing cos^ (HR) and sin^ (HR) in Eqs. and in terms of cos{2HR) and the 
subsequent Taylor expansion of cos{2HR) about R = would lead to the standard series 
expansion of the elements of the metric tensor in Fermi coordinates as in Eqs. ©^(H))- Such 
series are uniformly convergent for all R; however, the Fermi coordinates are admissible 
only for < 2HR < it. The hypersurface R = n/ {2H), where the timelike Killing vector Ot 
becomes null, is a static limit surface. 
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III. TIDAL DYNAMICS IN DE SITTER SPACETIME 

Let us now consider the radial motion of a test particle away from the reference observer 
O at -R = in the spherical Fermi coordinates We concentrate on the general reduced 
geodesic equation (jH)). In the case of radial motion in the spherically symmetric spacetime 
region of Eq. ()31|) , the only relevant nonzero connection coefficients are 

= -^Hsm{2HR), T^^ = T^^ = -Htan{HR). (32) 



Therefore, the radial equation of motion is 



Htan{HR) 



cos^HR) - 2 



dR'' ^ 
dT 



0. (33) 



This equation can be easily integrated once and the result for R = dR/dT is 

if = cos'^ (HR) -{1-V^) cos\HR), (34) 

where Vq = R{T = 0) is the initial speed at i? = 0. Next, reducing Eq. to quadrature 
and using the relation 

dx 1 , Vl + sin^ x + Vl + sin a; . . 



cosx sin^ X sin^ X — Vl + A^ sm x 

with A^ = — 1 + 1/Vq, it is possible to show that the general solution of the radial equation 
with the specified initial conditions is 

tan{HR) = ±Vosmh{HT). (36) 

In general, we require that the Fermi time T increase along the timelike worldline of every 
observer. Therefore, assuming Vq > 0, we are only interested in the upper sign in Eq. ()3(ij) . 
The modified Lorentz factors F and F^ are given in this case by 



To ^ To 



(37) 



cos^HR)' cos{HR)' 

where Fq = (1 — Vq)^^. The speed of the particle initially increases, remains constant 
or decreases depending upon whether Vq is less than, equal to or greater than the critical 
speed 1 / v^; nevertheless, Eq. (^Tj) indicates that regardless of the particle's initial speed, as 
T ^ oo it monotonically approaches a null geodesic at the static limit surface HR = 7t/2. 
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1.2 1.4 



FIG. 2: Plot of R versus HR, < HR < tt/2, for Vb = 0.4, 1/V2 and 0.9. In each panel, the lower 
(upper) curve represents the motion according to the exact (approximate) equation of motion in 
Fermi coordinates. In the exact case, R vanishes at HR = tt/2. 

Moreover, the speed of the particle, as "seen" by the reference observer, is zero at this 
boundary surface, where everything, including light, "appears frozen." We note that R has 
a maximum at cos{HR) = Fq/v^ for Vq < 1/v^ as illustrated in Figure El The case of a 
null ray follows from Eq. (jHU)) with the upper sign for Vq = 1; in fact, in this case Eq. ()36p 
can also be expressed as tanh(iyT) = sin{HR). 

Let us now determine to what extent these exact results are reflected in the first-order 
tidal terms of the generalized Jacobi Eq. (jllj) . For radial motion, we find 

d^R 



H^(1-2R^)R, 



(38) 
HR and 



which follows as well from Eq. (jHHI) . since to linear order we have ta.n{HR) 
cos{HR) ^ 1. Eq. (pH|l is thus valid for R sufficiently small compared to H^^. 

To study the general behavior of the solution of Eq. ()38j) . we note that this equation can 
be integrated once and the result is 
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where Vc = l/V^- Let us define a function z/) via tlie integral 



^^ix;u):= I -j=^=^, (40) 
'o V 1 + ve^^ 



where is a constant parameter. Near x = 0, we have the Taylor expansions 

*.(x;.)^(. + l)-V=|.T^|! + ^|! + ...]. (41) 

X = + ± + ^:<l!i±^4,| + ...]. (42) 

The general solution of Eq. (jHUjl is thus given by 

^4V2HR; 2Vo - 1) = ±HT. (43) 

For the physical problem under consideration here only the upper sign in Eq. (j43|) is needed; 
in this case, the motion can be described in terms of an attractor involving uniform motion 
at speed Vc. Explicitly, it is simple to see via the effective potential in Eq. fl39|) that for 
Vq > Vc, the particle monotonically decelerates and asymptotically (i.e. for HT —>■ oo and 
HR — > oo) approaches the critical speed. For Vq = Vc, the particle moves uniformly, while 
for Vq < Vc the particle monotonically accelerates and asymptotically approaches Vc. These 
results, together with the extent of the (initial) agreement between Eq. (j43|) and the exact 
result given by Eq. (j36|) . are illustrated in Figure |21 Analytic estimates for the difference 
between the exact solution and the approximation at a given time T can be obtained using 
Eqs. dnSl) and 



IV. FERMI COORDINATES IN GODEL SPACETIME 



Let us next consider the stationary and spatially homogeneous rotating universe model 



discovered by Godel 
can be expressed as 



17[. It has been discussed by a number of authors 



1811. Godel's metric 



ds"^ = -df - 2V2U{x)dtdy + dx^ - U^dy^ + dz^ , (44) 

where 

U{x) = e^^" (45) 
and n is a positive constant. The Ricci curvature for this spacetime is given by 

R^, = 2n\^u,, (46) 
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where = (5^o is the four-velocity vector field for free particles that are at rest in space 
and coincides with the timelike Killing vector field dt- The source of the Godel gravitational 
field could be a perfect fluid with velocity and constant density /i and pressure p given 
by = p = Vt^ / (Svr), where Vtdz is the vorticity vector of the geodesic worldlines of the fluid 
source. Alternatively, the source of the Godel fleld could be dust of constant density Vt^ / (47r) 
together with a cosmological constant A = — The investigation of the physical 
aspects of the Godel universe has led to the introduction of various interesting coordinate 
systems 0, |2l|; for instance, Gaussian coordinate systems for the Godel spacetime have 
been constructed in Ref. j^]. 

As in section InJ we are interested in the class of reference observers that are flxed in 
space; that is, they have constant x, y and z coordinates. These free observers are endowed 
with an orthonormal tetrad frame A'^j.^^ that is parallel transported along their geodesic 
worldlines. Thus A^^^^ = -u^ and A^^.^, 2 = 1,2, 3, can be expressed as 

^^(^i) = ^'^(^i) cos Qt + A'^j.g) sin fit, (47) 
A''^^) = sin nt + A''(2) cos (48) 

where in (t, x, y, z) coordinates we have ^^(^i-^ = S\ and 

A'^(,) = (-v^,0,^,o); (49) 

moreover, A'^^g-j = coincides with the spacelike Killing vector fleld d^- Thus at each 
point in space the unit gyro axes of the corresponding reference observer rotate about the 
z direction with frequency Q. For the sake of simplicity and without any loss in generality, 
we choose the reference observer atx = y = z = 0. 

It can be shown that all of the nonzero components of the Riemann tensor for metric (|44|1 
can be obtained from 

i?oioi = ^^', Ro202 = ^^U^ Rmi2 = -V2ffU, i?i2i2 = 3^2[/^ (50) 

via the symmetries of the Riemann tensor. It is then straightforward to demonstrate that 
the nonzero components of the projection of the Riemann tensor on the reference tetrad 
fleld can be obtained from 



(51) 
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via the symmetries of the Riemann tensor. Using Eqs. Q^Q, we note that in Fermi 
coordinates the spacetime metric — containing only the lowest-order tidal terms — would be 

ds^ ^-[1 + n^{X^ + Y^)]dT^ + dX^ + dY^ + dZ^ 

1 (52) 

- -n^xdv -Ydx)'^. 

o 

To find the exact expression for this metric, we need to specify all of the spacelike geodesies 
of the Godel spacetime that are orthogonal to the worldline of the chosen reference observer 
O : (t, x) = (t, 0). The spacelike geodesies of the Godel spacetime are given by 

t' + V2Uy' = E, V2Ut' + U'^y = k, (53) 

z' = h, -t'^ - 2V2Ut'y' + x'^ - U^y'^ + z'^ = 1, (54) 

where t' = dt/da, etc. Here a is the proper length of the spacelike geodesic and E,k and 
h are constants of integration. We need a general solution of these equations that would 
correspond to the geodesic segment from Q to P in Fig. 1. The vector tangent to this 
segment, namely, {f , x' , y' , z') at a = is at event Q : (t, 0); moreover, is orthogonal 
to A^^o), i.e. 

-^/^^(o) =t' + ^^Uy' = E = 0, (55) 

where Eq. fjS^ has been used. Therefore, with E = and z = ha, Eqs. fj3Hj) - (|M|) imply 
that 

t' = V2^, x'V^ = l-/.^ y' = -A (56) 
where h"^ < 1. The equation for x can be written as 

U'^ - 2(1 - h^)Q^U^ = -2Q^e, (57) 

which has the general solution 

U = a cosh(acr + b), (58) 



where a := Hy'2{l — h"^). Here a > and h are constants that are related by the requirement 
that t/ = 1 at cr = 0, i.e. 

a cosh 6 = 1. (59) 
Substitution of Eq. in Eq. (fFFj) results in 



aa 



V2n\k\. (60) 



14 



It follows from the solution of the other equations in ()56|) that 

2k 



t-T = (arctan e""^^ - arctan e*) , (61) 
k 

[tanh(a(T + 6) - tanh h] . (62) 



V2aQ\k 

Next, we compute '^^-^^(i) where a = 0, and use Eq. (0) to find Fermi coordinates 

(T, X, r, Z) such that 

T = T, X cos ^7T — y sin fiT = cr I /c I sinh 6, (63) 
X sin + F cos QT = -ak, Z = ah. (64) 

We note that + Y"^ + Z"^ = a^, since h"^ + k^/a'^ = 1. It proves useful to introduce 
cylindrical Fermi coordinates (p, if, Z) such that 

X = pcosip, Y = psinip. (65) 

Then, Eqs. can be written as 

k 

V2Qp = aa, cos((/9 + QT) = tanh 6, sin((/9 + QT) = ——a. (66) 

Consider now Eq. for x : expanding its right-hand side and using Eq. (jHUI) together 
with Eq. dSni), we find 

^V2nx ^ cosh(\/2(]p) + smh{V2np) cos((^ + QT). (67) 

Similarly, Eq. ()62|1 for ?/ can be written as 

^ tanh(V2np)sm(^ + nr) 

1 + tanh(y2fip) cos((^ + l^T) 



while Eq. flftlj) for t — T takes the form 



tan 



Ifi(T-t) 



1 - cos(^ + QT) + [1 + cos((/3 + (]T)]e^f^/'' 
It is advantageous to introduce new variables u and v by 



(69) 



u 



V2np, v = ^ + nT, (70) 
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so that the transformation to Fermi coordinates, {t, x, y, z) ^ (T, X, Y, Z), is given by 



tan 



-Q{T - 1) 



1) sinf 



1 — cos V + {1 + cos w)e" ' 
sinh u sin v 



V2ny 



z = Z. 



(71) 
(72) 
(73) 



cosh u + sinh u cos v ' 

It should be noted that transformations (1721) and (1731) bear a strong resemblance to those 
used by Godel to show the rotational symmetry of his metric about the z axis (see Ref. LJ, 
p. 449). 



One can show that 



dt + \Pl\J dy = dT + —(cosh u — 1) dv, 



dx^ + dy' 



2fi2 



{dv? + sinh^ u dv^) , 



so that the Godel metric 



ds"^ = -{dt + V2Udyf + dx^ + U^dy'^ + dz^ 



now takes the form 
ds^ = - 



+ — (cosh u — 1) dv 



+ 7^{dv'^ + sinh^ v dv'^) + dZ'^ 



(74) 
(75) 

(76) 

(77) 



in (T, M, f , Z) coordinates. From pdp = X dX + F dY and dip = X dY — F (iX, the form 
of the metric in Fermi coordinates (T, X, F, Z) is 



ds"^ = -(1 + L) c/T^ - 2nFdT{X dY -Y dX) + rfX^ + dY"^ + dZ^ 

+ -^{XdY -Y dXf, 
P 



(78) 



where p = a/X^ + 



^ / cosli w 1 

L{u) = - sinh^ m, -F(m) = ( 



u 



n{v) 



sinhw 



u 



l-2F{v). (79) 



The metric functions L, F and Ti depend upon coshw and cosh(2M); expanding these in 
uniformly convergent power series in -u = \/2VLp for any p > 0, we find the complete tidal 
expansion of the metric in Fermi coordinates. This result agrees with Eq. (jK^ to order Vt^. 

To examine the admissibility of Fermi coordinates in the sense of Lichnerowicz [4], we 
need to determine both {g^y) and {g^'^) in Fermi coordinates. Indeed, given our convention 
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that the signature of the metric is +2, Lichnerowicz admissibihty requires that the principal 
minors of these matrices be negative. We recall that for a symmetric n x n matrix M, the 
principal minors are defined by 



det 



(80) 



Mki ■■■ Mkk 

for = 1, . . . , ra. Employing (T, t>, p, Z) coordinates, it is possible to show that the admissi- 
bility conditions for the inverse metric are satisfied provided 1 + 7i, 

1 



i + n 



/ cosh-u 



(3 — coshw) 



(81) 



is positive. That is, the admissible Fermi coordinates must remain within a circular cylinder 
about the Z axis such that coshu < 3 or p < pmax, where 

V2, r:. 1.25 



Pmax = ln(l + 72) ^ — 



(82) 



Expressing the Fermi metric (j78j) in terms of cylindrical coordinates (T, p, y?, Z), we find 

ds^ = -(l + L) dT^ - 2Vtp^F dTdif + dp'^ + p'^{l + n) d^'^ + c/Zl (83) 

In this form, the Godel spacetime region under consideration is stationary and cylindrically 
symmetric, and the corresponding Killing vectors are Ot, d^p and dz- It is important to 
note that the square of the magnitude of the azimuthal Killing vector is p^(l + 7^), which is 
positive only for p < pmax- Thus for a given T and Z, the circle of radius p is spacelike within 
the admissible region, null at p = pmax and timelike for p > pmax- The admissible region 
therefore excludes closed timelike lines that are known to exist in Godel spacetime 

HQ. 

It is interesting to employ Fermi coordinates to study some of the properties of Godel 
spacetime. We recall that the source of the Godel solution has four-velocity that coincides 
with the timelike Killing vector in = [t, x, y, z) coordinates of Eqs. ()44|) - P^ . Under the 
transformation x^ ^ x'^ = (T, p, ip, Z), Uf^ = (9x"/9x^)mq,, so that using Ua = Qao we 

have 



dt 



u, 



V2U{ 



X] 



dy 



which can be simply computed using Eq. ()74|) . For the metric ^^i, given by Eq 
\/—g = {\/2il)^^ sinhu, and the only nonzero components of g^"^ are g^^ = g^^ = 1 and 



(84) 
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-TT 



X 



X + 



X 



X 



x' + i 



(85) 
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where x = coshw. It follows that = (1, 0, —Q, 0) in x'^ coordinates or, expressed geomet- 
rically, dt = dx — ^d^. 



The Gode 



nal paper 



spacetime is of Petrov type D and has five Killing vector fields. In his origi- 



Godel already discussed the simple symmetries associated with four Killing 
vectors, which in terms of {t,x,y,z) coordinates of Eqs. and are given by dt, dy, 
dz and — V2^lydy. The fifth Killing vector has a more complicated form 

K = 2V2U'^dt - 2V2nyd:, + {2^fy'^ - U-'^)dy. (86) 

This can be derived directly from the Killing equation; alternatively, one can use the simple 
Killing vectors in the new forms of the metric given in Eqs. ()77|) . (fTSj) or (jH^ to find it. For 
instance, one can use Eqs. (ffT|l - (f73|l to determine the partial derivatives of t, x, y and z 
with respect to v. The result is 

dv = ^iU-' - m - yd. + ^^[1 + {2nV - U^')]dy, (87) 

so that the "new" Killing vector is expressible in terms of those in the original coordinate 
system as 

2V2 nd^ = -2V2 dt + dy + K. (88) 
It is remarkable that pursuing such an approach to the symmetries of Godel spacetime. 



Obukhov already discovered the coordinate transformations fl71|) - (f7Sj) in Ref. 21 1. 



V. TIDAL DYNAMICS IN GODEL SPACETIME 

Let us now consider the motion of a free test particle in the exact Fermi coordinate 
system that we have constructed in Sec. IIVI Our purpose here is to compare the dynamics 
in the exact case with that using only the first-order tidal terms; that is, the solution of the 
equations of motion based on the metric form (j52|) . 

To simplify matters, let us limit our considerations to "radial" motion orthogonal to the 
Z direction. Specifically, we assume that at T = 0, (X, Y,Z) = and V = (Vq, 0, 0), where 
Vq G [0, 1). Starting from the exact Fermi metric (jH^ in cylindrical coordinates, we find the 
equations of motion 

dT 3-x' ^ ' 
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where, as in Sec. IIV| x = cosh.(\/2flp) and Tq is the Lorentz factor corresponding to Vq. Let 
w := X + I5 so that 



w = 2cosh^(^). (91) 
v2 

For < p < pmax, we have that 2 < w < 4. Equation (jH^ can be written in terms of w and 
reduced to quadrature as 

iii^^=== dx = ±V2 QT. (92) 



i: 



Next, formulas 2.261 on p. 81 and 2.266 on p. 84 of Ref. 
the general solution as 



can be employed to express 



1 vr 1 

-V2 arcsin Wi H , arcsin W2 + -{V2 ^=) = ±V2nT, (93) 

V2^ 2' ^J2^' 

where W\ and are given by 

V„=(x + 1) • 

It follows from Eq. (|89|) that in this "radial" motion with zero orbital angular momentum 
about the Z axis, p increases monotonically from zero and after reaching a maximum at 

returns to the origin in a time-symmetric fashion. For Vq : ^ 1, p : ^ Pmax! that is, 
a test particle remains within the Fermi frame, while a null ray reaches the null circular 
boundary. An interesting feature of p, derived from Eq. (|H^. should be noted: it vanishes 
at p = for Vo = 1/ v^, which is the critical speed for relative gravitational motion. 

The same kind of motion in the case of first-order tidal terms, metric ()52|) . turns out to 
be purely radial with the equation of motion 

p = -fi2(l-2p2)p (97) 

and initial conditions that at T = 0, p = and p = Vq. This equation is expected to be 
valid for Vtp sufficiently small compared to unity. Equation (jHTj) can be integrated once and 
the result is 

p2 = y^2_^^^2_^^2)g2nV. (98) 
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FIG. 3: Plot of p versus Vtp, Q<Vlp< ^/21n (1 + ^2), for Vq = 0.5, l/\/2 and 0.9. In the top and 
middle panels, the lower (upper) curves represent the motion according to the exact (approximate) 
equation of motion in Fermi coordinates. In the top panel, p monotonically increases in the 
approximate case. In the bottom panel, the inner curve corresponds to the exact equation of 
motion. 

moreover, the solution of Eq. ()98j) can be expressed as 

^+{V2Vtp- 2V^ - 1) = ±nT. (99) 

It is simple to see from the effective potential in Eq. (f^Hjl that for Vq > Vc, the particle 
monotonically accelerates to almost the local speed of light, while for Vq = Vc, the motion is 
uniform. For Vq < Vc, the motion along any radial axis in the {X, Y) plane is periodic and 
confined to the interval [—po,po], where 

We note that for < Vq < Vc, po > p, while for Vb = 0, po = P = 0. These results, together 
with a comparison with the exact solution, are illustrated in Figure El Analytic estimates 
for the difference between the exact solution and the approximation can be obtained using 
Eqs. and 
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VI. DISCUSSION 



The metrics of de Sitter and Godel spacetimes in Fermi coordinates can be expressed 
in infinite series of tidal terms that are uniformly convergent over all of space and time; 
however, the requirement of (Lichnerowicz) admissibility limits their respective domains of 
applicability to < HR < 7r/2 and < ^7p < \/2\n{l + \/2). In terms of the motion of free 
test particles in these domains, it turns out that — for the cases considered in this paper — the 
first-order tidal terms already provide good approximations over significant neighborhoods 
about the origins of these domains as illustrated in Figures 2 and 3. It is possible to provide 
detailed analytic estimates of the difference between the geodesic (deviation) equation and 
the generalized Jacobi equation using the solutions of these equations of relative motion 
presented in Sees. Ill and V. 



APPENDIX A: CRITICAL SPEED 



Imagine the radial motion of a swarm of particles away from a collapsed configuration as 
in an astrophysical jet. Relative to a free test observer and in a Fermi coordinate system 
based on this observer's worldline, free particles starting from the position of the observer 
and moving ultrarelativistically outward with speed above V"c = 1/ \/2 decelerate toward the 
critical speed l/\/2. However, free particles moving ultrarelativistically normal to the jet 
direction accelerate relative to the observer; via collisions with neighboring particles, the 
corresponding tidal energy can be imparted to the collapsed object's environment. One 
can in fact envision at the position of the observer a critical velocity cone with its axis 
along the jet direction and a total opening angle of 29c at its vertex such that tan 6'c = -\/2. 
Within this cone free ultrarelativistic particles decelerate relative to the observer, while 
they accelerate outside the cone. For infrarelativistic motion with speed below Vc = 1/ v^, 
tidal effects tend to behave more or less as one would generally expect on the basis of 
Newtonian gravitation theory. That ultrarelativistic gravitational tidal effects could exhibit 
astrophysically interesting features that would be contrary to Newtonian expectations was 
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JJ. It 



fi.t pointed out . aef. Q and Kas been tbe subject of seve.a. .eeent pape. 

is important to emphasize the general significance of the critical speed for relative motion 
in accelerated systems and gravitational fields. These results follow from the generalized 
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Jacobi equation that is based on the lowest-order tidal terms; we explain in Appendix B 
why it is impractical to use exact tidal terms for these astrophysically significant problems. 

It is interesting to point out here that in the specific context of (essentially radial) 
geodesic motion in the exterior Schwarzschild field, the critical speed was already discussed 
by Hilbert [23]. An examination of the early references to this subject is contained in 
Ref. j2J]. In terms of Schwarzschild coordinates, the critical speed is given by Vr = as 



26| derived the 



discussed in detail in Ref. |25[. Using a more invariant approach, McVittie 
critical speed K = 1/ a/2 for geodesic motion in Schwarzschild spacetime. 

Shapiro's observation of the gravitational time delay can be interpreted to imply that 
light slows down in the gravitational field of a massive body. This is in apparent confiict 
with the fact that in Newtonian gravity particles speed up as they fall toward a massive 
object. These ideas are properly integrated in general relativity through the concept of 



critical speed. In this context, Vc and Vc have been treated in Ref. [25|, which should be 
consulted for a more detailed treatment of (and further references to) this topic. Finally, we 
note that the critical speed Vc = 1/VS has been discussed recently in connection with the 
deflection of particles by a radially moving gravitational lens 



APPENDIX B: FERMI COORDINATES IN SCHWARZSCHILD SPACETIME 

The purpose of this appendix is to bring out the difficulties encountered in attempts to 
employ explicit Fermi coordinates in black-hole spacetimes. To this end, we imagine in the 
following treatment the simplest situation of astrophysical interest, namely, purely radial 
motion in the exterior Schwarzschild gravitational field. 

Consider the exterior Schwarzschild spacetime of a spherical source of mass M represented 
by the metric 

ds^ = _ (^1 _ ^) dt^ + (^1 _ ^ j dr^ + r\de^ + sin^ 9 d(j)^) (Bl) 

for r > 2M. To simplify matters, we limit our treatment to the motion of test particles 
along a fixed radial direction, which we can choose to be the z axis (i.e. 6 = 0) with no loss 
of generality. 

To construct Fermi coordinates for this two-dimensional world, we must specify a refer- 
ence observer. Imagine, therefore, a free observer O : {i, f) that starts from tq > 2M at 
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t = and follows an escape trajectory that reaches radial infinity with zero speed. The 
geodesic equations of motion for this observer are 



dt 1 df 2M 



(B2) 



r 

where r is the observer's proper time such that r = at t = and f = r^. The system flB2|) 
can be integrated and the result is 

f = 2M cosh^ w, St-AM cosh^ w = ci, (B3) 
3t - 4M cosh'^ w - 12M(cosh w + ln tanh - ) = C2, (B4) 

where Ci and C2 are constants of integration and can be expressed in terms of vq using the 
initial conditions. 

The observer carries an orthonormal parallel-propagated tetrad ^^(^a) ^-loiig its path. This 
local frame is given by unit vectors along the temporal and radial directions 



A%-((l-f)-V^.O,0). (B5) 



n3) = (\/^(l-^)-M.0.0). (B6) 

in {t, r, 6, 0) coordinates, respectively, as well as A'^j-j^-j and A'^^g) • The axial symmetry of the 
spacetime about the z axis implies that a rotational degeneracy exists in the choice of A'^j.^^^ 
and ^'^(^2)- -^^t' oiice these unit vectors are chosen at t = 0, they are then parallel transported 
along the path of O. 

Let us next consider the spacelike geodesies appropriate to the two-dimensional world of 
radial motion under consideration. They are given by 



dt p dr , 2M 

r 

where p is an integration constant and e = ±1. As in Figure ^ we need only those spacelike 
geodesies that are orthogonal to the reference worldline at Q. It follows from ^^A'^^qj = 
that 

Moreover, it follows from Eq. (Q) that in this case the Fermi coordinates are 

T = r, X = Y = 0, Z = ea. (B9) 



23 



It remains to integrate system ()B7|) from Q to P for constant T. This system can be written 

as 

dt p dr I 2 2M 

r 

where p = ep and g are functions of T, 



/2M / 2M , , 

p = \/— ^ = v^ + ~- ^^^^^ 



The solution of system ()B10|) is simphfied if we introduce a new quantity ( such that 



2M g 
Then, we find 



^coshC- (B12) 



C + I sinh(2C) -C-l sinh(2C) = ^Z, (B13) 
•^(C) - -^(0 + j(^-0 = ^[t- m - pZi (B14) 



where C, is given by 

,(1 

c, = cosn I 

and J^{C) is defined to be 



C = cosh-V^) (B15) 
p 



HO=H^^)- (B16) 
Here A and B are functions of T given by 

A = 1 + 2p2 + 2pg, 5 = 1 + 2p2 _ 2pg. (BIT) 

Inspection of Eq. ()B13j) makes it evident that it is not possible to express C, explicitly in 
terms of T and Z. It follows that the radial coordinate r in the Schwarzschild metric ()B1|) 
cannot be expressed explicitly in terms of T and Z. For the explicit expression of the 
metric in Fermi coordinates, one must resort to expansions in powers of Z. It is therefore 
clear that the transformation (t,r) i— > (T, Z) is given only implicitly by Eqs. ()B12|) - ()B17|) . 
Taylor expansions can be used to obtain an explicit form for the transformation to Fermi 
coordinates. We find that 

t = i{T) + hZ + 1/2^2 + y^Z' + 0{Z'), (B18) 
r = f(T) + Z + \^Z' - \^Z' + 0{Z% (B19) 
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where 



/i 



P 



f2 = - 



1 



h 



1 p^(4 — p^+p"' 



(B20) 



Using transformations ()B18|) - ()B19|) . we find the relevant metric coefficients in Fermi coor- 
dinates 

2M 



QTT = -1 + -^Z' + 0{Z'), qtz = 9ZT = 0{Z'), gzz = l + 0{Z'), 



in agreement with Eqs. (ID)-®, since in this case 

2M 



R 



■TZTZ 



(B21) 



(B22) 



We have shown that it is not possible to provide explicit Fermi coordinates for radial 
motion in the exterior Schwarzschild spacetime. The upshot of our direct approach is just a 
different approximation scheme for determining the series of tidal terms in Eqs. 0^(01) • In 
this way, we have alsoprovided the justification for using the generalized Jacobi equation 



in our previous work 
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APPENDIX C: EXTENT OF VALIDITY OF THE GENERALIZED JACOBI 

APPROXIMATION 

The generalized Jacobi equation (fTT|) is a first-order approximation to the geodesic devi- 
ation equation (jH} in Fermi normal coordinates. The purpose of this appendix is to describe 
a general method for determining the timescale over which the first-order approximation is 
valid. 

It is convenient to express Eqs. (jSj) and (jll|) in dimensionless form; to this end, we 
introduce T := T /TZ and A* := , where 7^ is a constant effective radius of curvature of 

the spacetime under consideration such that | A"! < 1 in the domain of admissibility of Fermi 
coordinates. In terms ofU := (A", V), Eqs. (jH)) and (fTT|) can be expressed respectively as 

W = (V,P(W) + Q(W)), (CI) 

U' = (V,P(W)), (C2) 

where a prime denotes differentiation with respect to T. Here 

p\x,v) := -7^2fi^o,o, + 2^i^.,,o^' 

+ li3^R..,V^V^ + ^R..,V^V^ + %..K^V^V')]^^ (C3) 
Q'{X,y) := ~[n{Tl^ -TlpV')V'^V^ + P\X,V)], (C4) 
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where V := dX'^/dT = (1, V). 

To estimate the difference between the solutions of Eqs. (jH)) and (fTT|). we introduce Z : = 
lA —U and note that 

Z' = J{U) - J{U) + SiU), (C5) 

where 

J(W) := (V, P(W)), S{U) := (0, Q(W)). (C6) 

The initial conditions for Eqs. (jH)) and (fTT|) are the same, hence it follows from the 
integration of Eq. ()C5|) that 

\ZiT)\ < r \J{U) - J{U)\ dT + r \S{U)\dT. (C7) 
Jo Jo 



It is clear from definitions ()C3|) and ()C4|) that Q consists of tidal terms of second, third and 
higher orders. Since \X\ < 1 and |V| < 1, we may suppose that ||-DJ|| < Aq and H^ll < Bq 
over an interval of interest [0, T] for some positive constants Aq and Bq. Hence it follows 
from Eq. that 

\Z{T)\ <Ao [ \Z{T')\ dT + BoT. (C8) 

n 

By an application of Gronwall's inequality [2^, we have the fundamental estimate 

|Z(T)| <fioTe^»^, (C9) 

which implies that 

\X{T) - X{T)\ < fioTe^"^. (CIO) 

As a measure of the validity of the approximation of A'(T) by X(T), we can determine 
the timescale % over which the magnitude of their difference is less than e, < e -C 1; that 
is, 

BoZe^"'^' = e. (Cll) 
Using Lambert's W function, which is the inverse of the function x ^ xe^, we find 



(C12) 



It should be clear that this method provides a general but crude estimate; indeed, sharper 
estimates may be obtained in specific cases using explicit solutions as in Sees. IIIII and El of 
this paper. 
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